Overview

The goal of this markdown is to see where in the Klamath basin, and what type of data is available from the WQX portal. The available databases on WQX are NWIS (USGS) and WQX(EPA). WQX is a data warehouse for water quality, biological, and physical data used by state environmental agencies, the EPA, other federal agencies, universities, private citizens, and others.

We are using two dataRetrival functions, whatWQPsites and readWQPdata

Detailed information about fields metadata can be found here

whatWQPsites

The following function returns sites that have collected data in the Klamath Basin. Below, is an overview of the column names of the data request are shown.

  • column names of site/station information
##  [1] "OrganizationIdentifier"                         
##  [2] "OrganizationFormalName"                         
##  [3] "MonitoringLocationIdentifier"                   
##  [4] "MonitoringLocationName"                         
##  [5] "MonitoringLocationTypeName"                     
##  [6] "MonitoringLocationDescriptionText"              
##  [7] "HUCEightDigitCode"                              
##  [8] "DrainageAreaMeasure.MeasureValue"               
##  [9] "DrainageAreaMeasure.MeasureUnitCode"            
## [10] "ContributingDrainageAreaMeasure.MeasureValue"   
## [11] "ContributingDrainageAreaMeasure.MeasureUnitCode"
## [12] "LatitudeMeasure"                                
## [13] "LongitudeMeasure"                               
## [14] "SourceMapScaleNumeric"                          
## [15] "HorizontalAccuracyMeasure.MeasureValue"         
## [16] "HorizontalAccuracyMeasure.MeasureUnitCode"      
## [17] "HorizontalCollectionMethodName"                 
## [18] "HorizontalCoordinateReferenceSystemDatumName"   
## [19] "VerticalMeasure.MeasureValue"                   
## [20] "VerticalMeasure.MeasureUnitCode"                
## [21] "VerticalAccuracyMeasure.MeasureValue"           
## [22] "VerticalAccuracyMeasure.MeasureUnitCode"        
## [23] "VerticalCollectionMethodName"                   
## [24] "VerticalCoordinateReferenceSystemDatumName"     
## [25] "CountryCode"                                    
## [26] "StateCode"                                      
## [27] "CountyCode"                                     
## [28] "AquiferName"                                    
## [29] "LocalAqfrName"                                  
## [30] "FormationTypeText"                              
## [31] "AquiferTypeName"                                
## [32] "ConstructionDateText"                           
## [33] "WellDepthMeasure.MeasureValue"                  
## [34] "WellDepthMeasure.MeasureUnitCode"               
## [35] "WellHoleDepthMeasure.MeasureValue"              
## [36] "WellHoleDepthMeasure.MeasureUnitCode"           
## [37] "ProviderName"

readWQPdata

Imports data from Water Quality Portal web service.

  • List of column names:
##  [1] "OrganizationIdentifier"                           
##  [2] "OrganizationFormalName"                           
##  [3] "ActivityIdentifier"                               
##  [4] "ActivityTypeCode"                                 
##  [5] "ActivityMediaName"                                
##  [6] "ActivityMediaSubdivisionName"                     
##  [7] "ActivityStartDate"                                
##  [8] "ActivityStartTime.Time"                           
##  [9] "ActivityStartTime.TimeZoneCode"                   
## [10] "ActivityEndDate"                                  
## [11] "ActivityEndTime.Time"                             
## [12] "ActivityEndTime.TimeZoneCode"                     
## [13] "ActivityDepthHeightMeasure.MeasureValue"          
## [14] "ActivityDepthHeightMeasure.MeasureUnitCode"       
## [15] "ActivityDepthAltitudeReferencePointText"          
## [16] "ActivityTopDepthHeightMeasure.MeasureValue"       
## [17] "ActivityTopDepthHeightMeasure.MeasureUnitCode"    
## [18] "ActivityBottomDepthHeightMeasure.MeasureValue"    
## [19] "ActivityBottomDepthHeightMeasure.MeasureUnitCode" 
## [20] "ProjectIdentifier"                                
## [21] "ActivityConductingOrganizationText"               
## [22] "MonitoringLocationIdentifier"                     
## [23] "ActivityCommentText"                              
## [24] "SampleAquifer"                                    
## [25] "HydrologicCondition"                              
## [26] "HydrologicEvent"                                  
## [27] "SampleCollectionMethod.MethodIdentifier"          
## [28] "SampleCollectionMethod.MethodIdentifierContext"   
## [29] "SampleCollectionMethod.MethodName"                
## [30] "SampleCollectionEquipmentName"                    
## [31] "ResultDetectionConditionText"                     
## [32] "CharacteristicName"                               
## [33] "ResultSampleFractionText"                         
## [34] "ResultMeasureValue"                               
## [35] "ResultMeasure.MeasureUnitCode"                    
## [36] "MeasureQualifierCode"                             
## [37] "ResultStatusIdentifier"                           
## [38] "StatisticalBaseCode"                              
## [39] "ResultValueTypeName"                              
## [40] "ResultWeightBasisText"                            
## [41] "ResultTimeBasisText"                              
## [42] "ResultTemperatureBasisText"                       
## [43] "ResultParticleSizeBasisText"                      
## [44] "PrecisionValue"                                   
## [45] "ResultCommentText"                                
## [46] "USGSPCode"                                        
## [47] "ResultDepthHeightMeasure.MeasureValue"            
## [48] "ResultDepthHeightMeasure.MeasureUnitCode"         
## [49] "ResultDepthAltitudeReferencePointText"            
## [50] "SubjectTaxonomicName"                             
## [51] "SampleTissueAnatomyName"                          
## [52] "ResultAnalyticalMethod.MethodIdentifier"          
## [53] "ResultAnalyticalMethod.MethodIdentifierContext"   
## [54] "ResultAnalyticalMethod.MethodName"                
## [55] "MethodDescriptionText"                            
## [56] "LaboratoryName"                                   
## [57] "AnalysisStartDate"                                
## [58] "ResultLaboratoryCommentText"                      
## [59] "DetectionQuantitationLimitTypeName"               
## [60] "DetectionQuantitationLimitMeasure.MeasureValue"   
## [61] "DetectionQuantitationLimitMeasure.MeasureUnitCode"
## [62] "PreparationStartDate"                             
## [63] "ProviderName"                                     
## [64] "timeZoneStart"                                    
## [65] "timeZoneEnd"                                      
## [66] "ActivityStartDateTime"                            
## [67] "ActivityEndDateTime"

Monitoring Sites

This table shows all the monitoring sites available across the Klamath Basin

Monitoring sites plot

Summary of sites by monitoring location type (MonitoringLocationTypeName: The descriptive name for a type of monitoring location)

Summary: Total number of records per organization

Summary: The top 50 Monitoring locations with the most data entries

Parameters

These are the water quality parameters available (pulling 2020 - 2024 data to view parameters). Parameters of interest will be explored on separate tabs

Temperature

Glimpse of temperature data pull

## Rows: 252
## Columns: 67
## $ organization_identifier                                <chr> "OREGONDEQ", "O…
## $ organization_formal_name                               <chr> "State of Orego…
## $ activity_identifier                                    <chr> "OREGONDEQ-2405…
## $ activity_type_code                                     <chr> "Field Msr/Obs"…
## $ activity_media_name                                    <chr> "Water", "Water…
## $ activity_media_subdivision_name                        <chr> "Surface Water"…
## $ activity_start_date                                    <date> 2024-05-07, 20…
## $ activity_start_time_time                               <chr> "14:44:00", "13…
## $ activity_start_time_time_zone_code                     <chr> "PST", "PST", "…
## $ activity_end_date                                      <date> 2024-05-07, 20…
## $ activity_end_time_time                                 <chr> "14:44:00", "13…
## $ activity_end_time_time_zone_code                       <chr> "PST", "PST", "…
## $ activity_depth_height_measure_measure_value            <dbl> NA, NA, NA, NA,…
## $ activity_depth_height_measure_measure_unit_code        <chr> NA, NA, NA, NA,…
## $ activity_depth_altitude_reference_point_text           <chr> NA, NA, NA, NA,…
## $ activity_top_depth_height_measure_measure_value        <dbl> NA, NA, NA, NA,…
## $ activity_top_depth_height_measure_measure_unit_code    <chr> NA, NA, NA, NA,…
## $ activity_bottom_depth_height_measure_measure_value     <dbl> NA, NA, NA, NA,…
## $ activity_bottom_depth_height_measure_measure_unit_code <chr> NA, NA, NA, NA,…
## $ project_identifier                                     <chr> "Surface Water …
## $ activity_conducting_organization_text                  <chr> NA, NA, NA, NA,…
## $ monitoring_location_identifier                         <chr> "OREGONDEQ-1076…
## $ activity_comment_text                                  <chr> NA, NA, NA, NA,…
## $ sample_aquifer                                         <chr> NA, NA, NA, NA,…
## $ hydrologic_condition                                   <chr> NA, NA, NA, NA,…
## $ hydrologic_event                                       <chr> NA, NA, NA, NA,…
## $ sample_collection_method_method_identifier             <chr> NA, NA, NA, NA,…
## $ sample_collection_method_method_identifier_context     <chr> NA, NA, NA, NA,…
## $ sample_collection_method_method_name                   <chr> NA, NA, NA, NA,…
## $ sample_collection_equipment_name                       <chr> NA, NA, NA, NA,…
## $ result_detection_condition_text                        <chr> NA, NA, NA, NA,…
## $ characteristic_name                                    <chr> "Temperature, w…
## $ result_sample_fraction_text                            <chr> NA, NA, NA, NA,…
## $ result_measure_value                                   <dbl> 11.1, 11.0, 19.…
## $ result_measure_measure_unit_code                       <chr> "deg C", "deg C…
## $ measure_qualifier_code                                 <chr> NA, NA, NA, NA,…
## $ result_status_identifier                               <chr> "Final", "Final…
## $ statistical_base_code                                  <chr> NA, NA, NA, NA,…
## $ result_value_type_name                                 <chr> "Actual", "Actu…
## $ result_weight_basis_text                               <chr> NA, NA, NA, NA,…
## $ result_time_basis_text                                 <chr> NA, NA, NA, NA,…
## $ result_temperature_basis_text                          <chr> NA, NA, NA, NA,…
## $ result_particle_size_basis_text                        <chr> NA, NA, NA, NA,…
## $ precision_value                                        <chr> NA, NA, NA, NA,…
## $ result_comment_text                                    <chr> NA, NA, NA, NA,…
## $ usgsp_code                                             <chr> NA, NA, NA, NA,…
## $ result_depth_height_measure_measure_value              <dbl> NA, NA, NA, NA,…
## $ result_depth_height_measure_measure_unit_code          <chr> NA, NA, NA, NA,…
## $ result_depth_altitude_reference_point_text             <chr> NA, NA, NA, NA,…
## $ subject_taxonomic_name                                 <chr> NA, NA, NA, NA,…
## $ sample_tissue_anatomy_name                             <chr> NA, NA, NA, NA,…
## $ result_analytical_method_method_identifier             <chr> "170.1", "170.1…
## $ result_analytical_method_method_identifier_context     <chr> "USEPA", "USEPA…
## $ result_analytical_method_method_name                   <chr> "Temperature", …
## $ method_description_text                                <chr> "https://www.ne…
## $ laboratory_name                                        <chr> "DEQ Laboratory…
## $ analysis_start_date                                    <date> 2024-05-07, 20…
## $ result_laboratory_comment_text                         <chr> NA, NA, NA, NA,…
## $ detection_quantitation_limit_type_name                 <chr> "Minimum Report…
## $ detection_quantitation_limit_measure_measure_value     <dbl> -10, -10, -10, …
## $ detection_quantitation_limit_measure_measure_unit_code <chr> "deg C", "deg C…
## $ preparation_start_date                                 <date> NA, NA, NA, NA…
## $ provider_name                                          <chr> "STORET", "STOR…
## $ time_zone_start                                        <dbl> 8, 8, 8, 8, 8, …
## $ time_zone_end                                          <dbl> 8, 8, 8, 8, 8, …
## $ activity_start_date_time                               <dttm> 2024-05-07 22:…
## $ activity_end_date_time                                 <dttm> 2024-05-07 22:…

Glimpse of station data pull

## Rows: 5,408
## Columns: 37
## $ organization_identifier                              <chr> "USGS-CA", "USGS-…
## $ organization_formal_name                             <chr> "USGS California …
## $ monitoring_location_identifier                       <chr> "USGS-11488495", …
## $ monitoring_location_name                             <chr> "LOST R NR HATFIE…
## $ monitoring_location_type_name                        <chr> "Stream", "Stream…
## $ monitoring_location_description_text                 <chr> NA, NA, NA, NA, N…
## $ huc_eight_digit_code                                 <chr> "18010204", "1801…
## $ drainage_area_measure_measure_value                  <dbl> NA, NA, 1.74, 9.9…
## $ drainage_area_measure_measure_unit_code              <chr> NA, NA, "sq mi", …
## $ contributing_drainage_area_measure_measure_value     <dbl> NA, NA, NA, NA, N…
## $ contributing_drainage_area_measure_measure_unit_code <chr> NA, NA, NA, NA, N…
## $ latitude_measure                                     <chr> "41.9538888900000…
## $ longitude_measure                                    <chr> "-121.50333330000…
## $ source_map_scale_numeric                             <chr> "24000", "24000",…
## $ horizontal_accuracy_measure_measure_value            <chr> "1", "1", "5", "5…
## $ horizontal_accuracy_measure_measure_unit_code        <chr> "seconds", "secon…
## $ horizontal_collection_method_name                    <chr> "Interpolated fro…
## $ horizontal_coordinate_reference_system_datum_name    <chr> "NAD83", "NAD83",…
## $ vertical_measure_measure_value                       <dbl> 4035.00, 4039.00,…
## $ vertical_measure_measure_unit_code                   <chr> "feet", "feet", N…
## $ vertical_accuracy_measure_measure_value              <dbl> 20.00, 20.00, NA,…
## $ vertical_accuracy_measure_measure_unit_code          <chr> "feet", "feet", N…
## $ vertical_collection_method_name                      <chr> "Interpolated fro…
## $ vertical_coordinate_reference_system_datum_name      <chr> "NGVD29", "NGVD29…
## $ country_code                                         <chr> "US", "US", "US",…
## $ state_code                                           <chr> "06", "06", "06",…
## $ county_code                                          <chr> "093", "093", "04…
## $ aquifer_name                                         <chr> NA, NA, NA, NA, N…
## $ local_aqfr_name                                      <chr> NA, NA, NA, NA, N…
## $ formation_type_text                                  <chr> NA, NA, NA, NA, N…
## $ aquifer_type_name                                    <chr> NA, NA, NA, NA, N…
## $ construction_date_text                               <chr> NA, NA, NA, NA, N…
## $ well_depth_measure_measure_value                     <dbl> NA, NA, NA, NA, N…
## $ well_depth_measure_measure_unit_code                 <chr> NA, NA, NA, NA, N…
## $ well_hole_depth_measure_measure_value                <dbl> NA, NA, NA, NA, N…
## $ well_hole_depth_measure_measure_unit_code            <chr> NA, NA, NA, NA, N…
## $ provider_name                                        <chr> "NWIS", "NWIS", "…

Looking how many data entries there are for each day in 2024

Plot below shows the same data than table. This indicates that data is not being collected every day

Temperature from 2020 to 2025

Map below shows where monitoring locations are, and what agency is monitoring/maintain it

## Joining with `by = join_by(organization_identifier, organization_formal_name,
## monitoring_location_identifier, provider_name)`

Table below shows the amount of data entries per site, per year and month

Table below shows the amount of data entries per site, per year (note that this is per each site)

Table plot below shows the amount of data entries per organization, per year (note that this is per each organization, which includes all monitoring stations managed by them)

## Warning: Returning more (or less) than 1 row per `summarise()` group was deprecated in
## dplyr 1.1.0.
## ℹ Please use `reframe()` instead.
## ℹ When switching from `summarise()` to `reframe()`, remember that `reframe()`
##   always returns an ungrouped data frame and adjust accordingly.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's 
##   -0.64   13.12   16.06   16.30   19.48   73.09       8

Dissolved oxygen (DO)

Glimpse of dissolved oxygen data pull

## Rows: 3,832
## Columns: 67
## $ organization_identifier                                <chr> "HVTEPA_WQX", "…
## $ organization_formal_name                               <chr> "Hoopa Valley T…
## $ activity_identifier                                    <chr> "HVTEPA_WQX-TR_…
## $ activity_type_code                                     <chr> "Field Msr/Obs-…
## $ activity_media_name                                    <chr> "Water", "Water…
## $ activity_media_subdivision_name                        <chr> NA, NA, NA, NA,…
## $ activity_start_date                                    <date> 2020-06-02, 20…
## $ activity_start_time_time                               <chr> "00:00:00", "00…
## $ activity_start_time_time_zone_code                     <chr> "PST", "PST", "…
## $ activity_end_date                                      <date> 2020-06-02, 20…
## $ activity_end_time_time                                 <chr> "23:00:00", "23…
## $ activity_end_time_time_zone_code                       <chr> "PST", "PST", "…
## $ activity_depth_height_measure_measure_value            <dbl> NA, NA, NA, NA,…
## $ activity_depth_height_measure_measure_unit_code        <chr> NA, NA, NA, NA,…
## $ activity_depth_altitude_reference_point_text           <chr> NA, NA, NA, NA,…
## $ activity_top_depth_height_measure_measure_value        <dbl> NA, NA, NA, NA,…
## $ activity_top_depth_height_measure_measure_unit_code    <chr> NA, NA, NA, NA,…
## $ activity_bottom_depth_height_measure_measure_value     <dbl> NA, NA, NA, NA,…
## $ activity_bottom_depth_height_measure_measure_unit_code <chr> NA, NA, NA, NA,…
## $ project_identifier                                     <chr> "EPACONTINUOUS;…
## $ activity_conducting_organization_text                  <chr> NA, NA, NA, NA,…
## $ monitoring_location_identifier                         <chr> "HVTEPA_WQX-TR_…
## $ activity_comment_text                                  <chr> NA, NA, NA, NA,…
## $ sample_aquifer                                         <chr> NA, NA, NA, NA,…
## $ hydrologic_condition                                   <chr> NA, NA, NA, NA,…
## $ hydrologic_event                                       <chr> NA, NA, NA, NA,…
## $ sample_collection_method_method_identifier             <chr> NA, NA, NA, NA,…
## $ sample_collection_method_method_identifier_context     <chr> NA, NA, NA, NA,…
## $ sample_collection_method_method_name                   <chr> NA, NA, NA, NA,…
## $ sample_collection_equipment_name                       <chr> NA, NA, NA, NA,…
## $ result_detection_condition_text                        <chr> NA, NA, NA, NA,…
## $ characteristic_name                                    <chr> "Dissolved oxyg…
## $ result_sample_fraction_text                            <chr> NA, NA, NA, NA,…
## $ result_measure_value                                   <chr> "9.64", "8.6433…
## $ result_measure_measure_unit_code                       <chr> "mg/L", "mg/L",…
## $ measure_qualifier_code                                 <chr> NA, NA, NA, NA,…
## $ result_status_identifier                               <chr> "Final", "Final…
## $ statistical_base_code                                  <chr> "Maximum", "Mea…
## $ result_value_type_name                                 <chr> "Calculated", "…
## $ result_weight_basis_text                               <chr> NA, NA, NA, NA,…
## $ result_time_basis_text                                 <chr> "24 Hours", "24…
## $ result_temperature_basis_text                          <chr> NA, NA, NA, NA,…
## $ result_particle_size_basis_text                        <chr> NA, NA, NA, NA,…
## $ precision_value                                        <chr> NA, NA, NA, NA,…
## $ result_comment_text                                    <chr> NA, NA, NA, NA,…
## $ usgsp_code                                             <chr> NA, NA, NA, NA,…
## $ result_depth_height_measure_measure_value              <dbl> NA, NA, NA, NA,…
## $ result_depth_height_measure_measure_unit_code          <chr> NA, NA, NA, NA,…
## $ result_depth_altitude_reference_point_text             <chr> NA, NA, NA, NA,…
## $ subject_taxonomic_name                                 <chr> NA, NA, NA, NA,…
## $ sample_tissue_anatomy_name                             <chr> NA, NA, NA, NA,…
## $ result_analytical_method_method_identifier             <chr> NA, NA, NA, NA,…
## $ result_analytical_method_method_identifier_context     <chr> NA, NA, NA, NA,…
## $ result_analytical_method_method_name                   <chr> NA, NA, NA, NA,…
## $ method_description_text                                <chr> NA, NA, NA, NA,…
## $ laboratory_name                                        <chr> NA, NA, NA, NA,…
## $ analysis_start_date                                    <date> NA, NA, NA, NA…
## $ result_laboratory_comment_text                         <chr> NA, NA, NA, NA,…
## $ detection_quantitation_limit_type_name                 <chr> NA, NA, NA, NA,…
## $ detection_quantitation_limit_measure_measure_value     <dbl> NA, NA, NA, NA,…
## $ detection_quantitation_limit_measure_measure_unit_code <chr> NA, NA, NA, NA,…
## $ preparation_start_date                                 <date> NA, NA, NA, NA…
## $ provider_name                                          <chr> "STORET", "STOR…
## $ time_zone_start                                        <dbl> 8, 8, 8, 8, 8, …
## $ time_zone_end                                          <dbl> 8, 8, 8, 8, 8, …
## $ activity_start_date_time                               <dttm> 2020-06-02 08:…
## $ activity_end_date_time                                 <dttm> 2020-06-03 07:…

Looking how many data entries there are for each day from 2020 to today

Plot below shows the same data than table. This indicates that data is not being collected every day

Map below shows where monitoring locations are, and what agency is monitoring/maintain it

Table below shows the amount of data entries per site, per year and month

Table below shows the amount of data entries per site, per year (note that this is per each site)

Table plot below shows the amount of data entries per organization, per year (note that this is per each organization, which includes all monitoring stations managed by them)

## Warning: Returning more (or less) than 1 row per `summarise()` group was deprecated in
## dplyr 1.1.0.
## ℹ Please use `reframe()` instead.
## ℹ When switching from `summarise()` to `reframe()`, remember that `reframe()`
##   always returns an ungrouped data frame and adjust accordingly.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's 
##   0.000   8.525   9.348  18.749  10.200 943.000      55